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This study investigates the modifications of the thermal spectrum, at low frequency, 
induced by an external damping on a system in heat contact with internal fluctuat- 
ing impurities. Those impurities can move among locations and their oscillations are 



associated with a loss function depending on the model. 



The fluctuation properties of the system are provided by a potential function shaped 



^ , by wells, in such a way that jumps between the stationary positions are allowed. The 

O 

" ' power spectral density associated with this dissipation mechanism shows a /~ 7 tail. 

I The interest of this problem is that many systems are characterized by a typical 

/~ 7 spectral tail at low frequency. The model presented in this article is based on a 
threshold type behaviour and its generality allows applications in several fields. 



The effects of an external force, introduced to produce damping, are studied by 

o ' 

using both analytical techniques and numerical simulations. 



The results obtained with the present model show that no reduction of the power 
spectral density is appreciable below the main peak of the spectral density. 
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I. INTRODUCTION 



Since gravitational wave interferometers will sense the passage of a gravitational wave as a difference 
in the light path between mirrors in the two arms, every kind of excitation is a noise source. Cooling the 
mirrors of these detectors, by measuring and controlling with a feedback loop the principal fundamental 
motions that may be thermally excited, may offer a means to reduce thermal noise. 

With this motivation in mind, we investigated models that generate frequency spectra with a tail 
for / — > 0. This dependence is found in many completely different physical systems, suggesting 
a possible underlying simple mechanism. The typical physical parameters are different for each of 
them, and especially the characteristic microscopic quantities. But many classes of phenomena are 
characterized by the same value of 7. 

Indeed the responsive behaviour of every system is classified according to the parameters of the 
phenomenological equations by which the main characteristics of its evolution and interactions are 
described. Moreover the contribution of the microscopic degrees of freedom is absorbed in the definition 
of the parameters so that it is possible to describe the motion of a system in terms of a phenomenologial 
equation for the macroscopic degrees of freedom. The microscopic degrees of freedom motions are 
integrated over as in averages and their effect in the dynamics of the system is contained in the physical 
constants. For example in high energy physics the masses are corrected for those particles which are 
not isolated but interacting with other ones, and since measurements imply interactions, the measured 
values are the renormalized values instead of the bare values. Similarly applications of those principles 
may be found in low energy physics fUJ^]. A typical example is the up-conversion of resonant frequencies, 
related with a non-linear perturbation term that affects both the bare value of the resonant frequency 
and its width ||. 

In this paper a very simple model is studied whose basic characteristic is the switching between 
equilibrium configurations across energy barriers. There are no other requirements and this generality 
is the main feature of this flexible model that can be easily adapted for a lot of different physical 
phenomena. The principal interest has been the investigation of the influence on the spectral function 
and particularly on the typical low frequency tail of a damping force. 
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II. THE MODEL 



For those materials characterized by a quality factor that may be increased by decreasing the 
amount of impurities, a reasonable picture is that the motion of those degrees of freedom may generate 
a fluctuation-dissipation process across the sites in the crystal. In other words the impurities are 
supposed to move following a random walk. 

The energy of those impurities is related with their temperature; in cryogenic conditions they hardly 
can pass through the potential barrier and change their location. The probability of a jump from one 
site to another depends upon T. The system is characterized by its statistical properties and the 
temperature is defined as a stationary thermodynamical quantity. 

A simple random walk model may be considered to represent the switching between the equilibrium 
positions for a two well potential function 

It has two minima for x = ±\f \ pointed out in Fig.|]. 

Increasing the constant a implies a smaller distance between the two wells. From the Taylor expan- 
sion the following expression is achieved 

V{x)=mul[-^ + {x±^ r f + ..] (1) 
4a yet 

for x near the two minima. The constant ttuOq has been introduced so that has the physically 

right dimension of frequency. In the phase space 

V 

x = — 
m 

p = mujQX^l — ax 3 ) 

whose solutions represent the trajectories that depend on the initial conditions. 
The problem may be treated perturbatively near the two minima. 

Using the linear approximation for small oscillations a stochastic force may be introduced that 
satisfies the fluctuation-dissipation theorem. 
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In other words, the deterministic part of this force and the random component are related in such 
a way that the exchanged energies are balanced. 
The resulting equation is 

mx + — °x + 2mu 2 (x ± -= = J t < WW) >= W ~ 

<y v« V Q 

near the two minima ^fa~^. The constant Q is so large that underdamped conditions are always 
satisfied. 

In order to take into account the ratio of the thermal equilibrium energy over the local maxima in 
the potential function, a scale factor may be used in order to recover the equation 



where x 1 




As a consequence, a, that had the same dimensionality as x 2 , has become dimensionless. 
It is now easy to see that flat wells and high temperature are equivalent and the constant a' = ct^jr 
sums up the two effects. 

All the comments referred to Fig.|l| are valid for a' as well. 

More intutively the scaling of the x-axis by squeezing or stretching gives the potential function a 
different shape corresponding to making the wells nearer or farther, deeper or flatter. 
Following an expansion of the equation of motion near ^fa'~^ 

x'tf) = G(f)[u ^af) T 3V^u 2 J^df'x'(f')x'(f - /') + 

/oc roc 
df / df"x'(f')x'(f")x'(f -f- /")] 
-00 J —oo 

where the origin of x' has been fixed in ±^7=7- It may be noted that x' appears also on the right side 
of the equation and this means that this definition is implicit. Using a perturbative expansion gives 

Af) = <U) + Gif^o^M T 3vW/^d/Y (/Vo(/ - /') + 
df / df"x' (f')x' (f")x' (f -f- /") + 

-OO J —OO 

/oo poo 
d/V (/ - f')u 2 G(f) / df"x' (f")x' (f - f") + ...} 
-00 J— 00 
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where x' Q represents the solution in the zero order approximation. 

If x' Q is considered and the two point correlation function is calculated, the associated integral over 
frequencies is 

/OO /'OO /'OO poo I 

df / df < x' (f)x' (f) >= df df'S(f)5(f + f') = - 
-OO J — oo J —oo J —oo £ 

for each of the two linearly approximated solutions near the minima. 

The function S(f) is called power spectral density. Taking into consideration higher order terms 
for x' gives corrections to S(f). For the property < £ >= the term representing the first correction 
of S(f) is proportional to a'. It is 

<x'(f)x'(f')>-5(f + f)S(f) = 

/oo 
d/i5(/i)5(/ - fi) + 
-oo 

/oo 
dfiS(h) + 
-oo 

/oo 
d/iS(/i) + 
-oo 

/oo 
d/i5(/i)»[G(/)G(/ - /i)] + . . .] 
-oo 

where G(f) = (-u 2 + ^ + 2wg)" 1 . 

All the terms may be easily derived using the Feynman diagram technique. In this paper indeed the 
analogies between a wave equation and a heat equation are used and the solution x'(f) is represented 
and expanded in a graphic way. The analogy with the Feynman diagram techniques is defined by 
associating a line with the Green's function G(f) and a cross with the force £(/). When two crosses are 
combined together a 5(f+f) arises. The iterative method of getting the terms of the expansion consists 
of substituting for the solution in the implicit definition at any order, the expression corresponding to 
the lower one in a'. 

When a' = the power spectral density may be calculated for each one of the two linear approxi- 
mations 

S(f) = CG(f)G(-f) C = ^ . 
Integrating over all frequencies S(f) gives \ for each of the two zero order expansions. This is due 
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to the fluctuation-dissipation theorem that links the dissipative force to the constant in front of £ in 
the motion equation. 

The Feynman diagrams are shown in FigfJ with the graphic representation of the solution. When 
two x'(t) are connected, the result is a tadpole type diagram. Indeed the mean value of x' 2 (t) does not 
depend on t if the conditions are stationary. In this context it should be reminded that the fluctuation- 
dissipation mechanism has been extended to the perturbed case. This generalization is to be intended 
as a limit of the methods and it is based on the physical assumption that the results are slightly modified 
as it is assumed in field theory techniques. 

The final result is 

< x (f)x U ) > =<>(/ + / J^; r~2 H 

Q (^ 2 - 2a; 2 ) 2 + ig# (to* - 2u, 2 ) 2 + ^ ~ 



(u;2 + g)[(^-8a; 2 )2 + ^ 



T2 



(- 2 + ^)(- 2 + *-64c 4 -^ 4 



to 2 - Sul + & 



+ 



' f >2 

+ 1^(7$ - 3^ 2 ) + a; 2 (10aA; 2 - 16u; 4 



3M/',.,2 _ ,.,2\ , ^o/^ 2 



(a; 2 - 2a; 2 ) 2 + ^ 

where the first line represents the graphs (0),(2) and (3) in FigfJ and the corresponding contribution 
modifies the definitions of the parameters but does not introduce a different dependence on time. 
This means that the typical shape of the power spectral density is not changed excepted some small 
corrections of its parameters. 

The characteristics of the spectrum are indeed changed by the other contributions involving con- 
volution integrals ||. These terms correspond to the graphs (1) and (4) that introduce a new feature 
in the shape of the spectrum. Indeed, they modify the power spectral density at low frequency for 
7^ << u) « ojq. In this range the spectrum decreases with increasing oj. The spectral density changes 
its shape when the value of various parameters is modified, but there is still a tail for that interval 
falling as o>~ 7 for < 7 < 2. 

Moreover a new feature accounted for in the formula above is a new peak that is produced at a 
frequency which is approximately twice the value of the main peak frequency. All these analytical 
results are obtained using a perturbative method that assumes that the system is never very far from 




the equilibrium conditions. 

If the system is allowed to go back and forth between the two minima of the potential function, an 
alternative approach is needed in order to simulate the dynamics of the system. A simulation code has 
been carried out and used to analyse different interesting cases. 

The physical picture is that the impurity has some chance of oscillating between two equilibrium 
positions acquiring energy from thermal excitations. 

It may be pushed into one direction or the other isotropically and the representation of those 
fluctuations is simply given by the stochastic term. 

If some kicks are in the same direction, energy is gained to overcome the potential barrier. This 
model reproduces the competition between thermal excitation and maxima between wells in the po- 
tential function. If the temperature is high the passage through the barrier has more probabilities to 
occur as it might be expected since T is a statistical definition of energy. 

At this point there is no straightforward way for keeping the damping term in the non-linear regime, 
apart from making physical guesses about the type of damping process. Choosing other forms would 
mean different spectral densities since the dynamics would be different. 

The relation between the random force and the damping term that is defined by the fluctuation- 
dissipation theorem in the linear approximation, has been used in the stochastic equation outside 
the two regions near the minima as well. It provides a solution whose distribution in stationary 
conditions satisfies the Boltzmann one. In general the Langevin equation implemented in the simulation 
code corresponds to a Fokker-Planck equation for the probability distribution that has the Boltzmann 
distribution as its stationary solution. This is one of the motivations to keeping the same expressions 
for the damping force and the random term in the non-linear case. The other motivation is that the 
nature of the damping effect is often defined by the properties of the heat-bath rather than on the 
mechanical characteristics of the physical system in contact with it. Another point worthy of noting 
is that the natural scale for the x-axis is related with the magnitude of the fluctuating forces. This is 



why the variable x' — x 




7^ has been used throughout our analyses. 
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III. SIMULATION AND NUMERICAL RESULTS 



The analogy existing between the Schroedinger equation and the Fokker-Planck equation can be 
exploited and an algorithm may be written based on this similarity. In particular the drift term and the 
diffusion term in the Fokker-Planck equation, play the same role of the potential function and kinetic 
term in the Hamiltonian operator. 

Using the same rules for path-integral calculations, the evolution operator may be factored in such 
a way that successive steps are performed which make the system fluctuate according to the diffusion 
term or make a step according to the deterministic term. 

Depending upon the required precision of the stepping in the simulated states of the system, a 
particular number of intermediate steps is needed to make the system evolve from its state to the 
successive one. This technique is simply an application of the Baker-Campbell-Hausdorff formula to a 
stochastic phenomenon @||. 

Recall that large values of ol imply high T or flat wells. In this high energy regime, the dynamics 
at low frequency are such that there is a constant trend as in harmonic oscillators damped by a viscous 
force. In Fig.|3| the spectral density is shown for a' = 10. No tail is appreciated at low frequency. 

If the value of a' is decreased the shape of the potential function is deeply characterized by two well 
defined minima and this affects the dynamics at low frequency where the trend of the power spectrum 
S(f) is not constant. On the contrary there is a slope at low frequency that in Fig.[| is quite visible 
and fitted with /~ 13 for the specific case of a' = 1. Indeed according to the tests I made using different 
values for the parameter a' the low frequency tail has the form J -7 with < 7 < 2. 

The dynamical model we constructed shows a low frequency tail in its spectrum. How is this system 
affected by a frequency independent damping force? Will the low frequency tail be attenuated by this 
external force? This system is a good candidate for investigation as the main fundamental modes are 
easily detected and suppressed by actively feeding back to the resonant motion. The damping we will 
use is an external viscous force proportional to x'. Our aim is to find out whether an external damping 
mechanism can suppress S(f) at low /. 

Some problems are to be faced in numerical studies. The numerical fluctuations in simulations are 
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the statistical errors that affect the mean value. They are reduced by averaging over many sets of data 
in order to reduce the deviations 

The statistical fluctuations are especially large at low frequency so that if we are interested in 
obtaining a smooth curve with small error bars, many runs must be averaged over. 

In Fig.|5] the spectral density is shown when the damping is applied. Only one of the peaks is well 
defined and indeed it becomes sharper since the damping force limits the motion to small oscillations 
around one of the minima of the potential function. 

The situation at low frequency is more interesting. The values at / ~ are not affected. In the 
range below the peak the low frequency tail is changed as /~ 2 . A similar trend may be obtained by 
increasing the depth of the wells so that the motion is mostly confined to oscillations inside one of the 
wells. 

Many models have been proposed which involve a tail f~ l , as it is typical of electronic circuits' 
noise. 

For mirrors the value of 7 for a spectrum /~ 7 is not certain. In fact, making experimental mea- 
surements of the fluctuations at low frequency is difficult even for state-of-the-art optics reasearch. 
Experimental data are mostly based upon resonance measurements and spectral curves are extrapo- 
lated for very low frequencies. 

For lack of off-resonance measurements linear equations have been studied with a complex elastic 
force 

The random-walk problem studied in this work focuses on only one property that is the existance 
of more than one equilibrium configuration. A similar model is used to reproduce the switching of 
polymeric units between two positions which are energetically equivalent. 

The simulation program has also been modified in order to investigate what may happen if there 
are more than two minima for the potential function. Instead of the two wells, a periodic potential has 
been studied giving the same results. In Fig.|6| the spectra obtained for two values of a' are shown. The 
fluctuation-dissipation theorem is extended to the general case although it could only be applied in the 
linear approximation near one of the minima of the potential function. Also x may be scaled in order 
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to have x' in the potential function 

V(x') = u>g(l -cosv^x') 

From the physical point of view a potential barrier is needed to represent the limits of the periodic 
structure that otherwise would be infinitely large, and this is not realistic. 

In this case a normalized stationary solution of the Fokker-Planck equation is the corresponding 
Boltzmann distribution which takes into account the modified potential function. In other words what 



is needed is a differentiable term that provides V(x') — > oo for both the limits x' — > ±00 |14 . 

We assumed that the size of the periodic structure is so large that no border condition affects the 
model. 



IV. CONCLUSIONS 

This unidimensional model is characterized by a tail at low frequency in its power spectral density 
that is not modified even when a frequency independent damping force is applied. 

The data are obtained from a time domain simulation. We have considered both a simple case for a 
two wells potential function and the more general case of a periodic frame with oscillations between the 
minima. The fluctuating motion of the degree of freedom and the associated dissipation are a model 
to produce a typical low frequency tail /~ 7 in the power spectral density S(f). 

The slope of such low frequency tail can be tuned by varying the constants; it depends on the 
competition between the temperature and the depth of the wells. 
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FIG. 1. The potential function ~ — ^- + a^-. The wells are farther apart if a decreases. Their depth depends 
on a as well, and scaling x so that every quantity becomes dimensionless gives a' ~ aT. Thus if T increases 
the results are equivalent to flattening the wells 
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FIG. 2. The perturbative corrections to the power spectral density. The graphs are obtained using the 



Feynman formalism for the correlation function, similar to a "propagator" 



Numerical results for a' = 10 and no additional external forces 
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FIG. 3. The power spectral density is shown for a' = 10. Many values have been tested for loq and Q which 
only influence the shape of the peaks. For this simulation the values chosen are u>o = 2.7Hz and Q = 100 
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FIG. 4. The power spectral density is shown for ol = 1, Uq = 1.1Hz and Q = 100 and these results can be 
compared to the ones in Fig|3|. The two wells are deeper and the curve gains a tail at low frequency that is 
/~ L3 . After reaching stationary conditions the averages over 5000 runs were computed 
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FIG. 5. An external damping force is introduced and its impact on the spectral curve is concentrated in the 
range of frequency around the peaks. Indeed one of them is substantially depressed 
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FIG. 6. Another example of tuning a constant to have a low frequency tail. The deterministic force is periodic 
on a lattice unidimensional space. The same scaling law applied previously, that takes the temperature into 
account, is used here to make quantities dimensionless 
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